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Résumée 


Caractérisation Thermo-optique d f un Distillateur Thermosolaire à 
Effet-Serre à Multiples-Effet sur site d'Adrar 

Résumée: 

La distillation des eaux saumâtres par effet-serre est limitée par la fraction de 
la radiation solaire incidente sur l'ouverture du capteur et par la chaleur latente 
d'évaporation étroitement liée au gradient de température entre la source chaude 
(la saumure) et la source froide (le condenseur). Le travail entrepris dans ce mémoire a 
permis de mieux comprendre les bilans énergétiques mis enjeu avec l'ambiance dans le 
distillateur à effet serre à multiple effets, et les propriétés thermodynamiques 
d'évaluation. La caractérisation thermo-optique du système est effectuée par une 
méthode simple basée sur des bilans énergétqiues et des propriétés thermo-optiques 
sous des conditions de fonctionnement constantes. La détermination du rendement 
optique et du facteur des pertes thermiques est effectué par la méthode des moindres 
carrées sur le nuage de points de mesure entre le rendement à la température réduite. 
Les résultats de caractérisation du distilalteur à multiple effets sur site d'Ardar montrent 
que le rendement optique et le facteur des pertes thermqiues globales obtenus sont 
(a,~80%, U l ~ 6.2 W/m 2 C) sous les conditions les plus proches du standard ASHRAE. 
Ces résultats sont relativement acceptables en comparaison avec les caractéristiques 
calculées qui sont (ctt~87%, Ul~ 5,5 W/m 2 C) et tenant compte des hypothèses 
simplificatrices du modèle et des ererurs de mesure. Le distillateur à multiples effets 
permet d’améliorer la production du distillât de 10 à 24% par rapport au distillateur à 

simple effet. 


Mos-clés: Distillation solaire, rendement optique, facteur des pertes thermiques, 
caractérisation, distillateur solaire à effet de serre multiples effets. 
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Introduction générale 


Introduction générale: 

La demande en eau douce en Algérie, croit chaque année de 4 à 5%, tandis que les 
ressources naturelles restent invariables. Le problème majeur, qui se pose aujourd'hui à 
l'homme avec une acuité de plus en plus accrue n'est certainement par d'ordre 
quantitatif, mais plutôt de nature qualitative. En effet, les eaux de qualité requise sont de 
plus en plus rares même dans les pays relativement bien dotés par la nature, suite à un 
processus de développement hautement polluant surtout dans le monde 
industrialisé[12]. La distillation solaire est un procédé de production d’eau douce à fort 
économique, surtout dans les régions caractérisées par un bon ensoleillement[3], 

Les régions sahariennes en général et la région d'Adrar en particulier recèlent 
d'importants gisement hydriques et solaires qui dépassent les 2200 kWh/m 2 /an sur une 
superficie aplanie de 428000 m 2 , très favorable pour la production de l’eau distillée à 
grande échelle permettant de réduire les frais de transport et de stockage et le 
développements des industries cosmétiques, chimiques à valeur ajoutée. 
Il est communément admis que la distillation des eaux saumâtres par effet-serre est 
limitée par la fraction de la radiation solaire incidente sur 1 ouverture du capteur et par la 
chaleur latente d'évaporation étroitement liée au gradient de température entre la source 
chaude et la source froide. En basse et moyenne température, le transfert de chaleur est 
une propriété de surface. Cependant, les distillateurs à multiple effets basés sur les 
propriétés optiques et thermiques permettent d'améliorer le rendement par rapport aux 
distillateurs à simple effet à effet serre. La quantité de distillât produite est fonction de 
la qualité de design du système et de l'énergie disponible sur site. Ainsi, la connaissance 
des performances thermiques et optiques de fonctionnement optimal du système sous 
des conditions de fonctionnement standards est cruciale pour décrire, comprendre et 
optimiser son design. 

L'objectif de ce travail intitulé: Caractérisation Thermo-Optique d'un Distillateur 
Thermosolaire à Effet-Serre à multiples-effet sur site d'Adrar, est la détermination de 
rendement optique et le facteur des pertes thermiques de système. On décomposé ce 
travail en troix chapitre: 

Le premier chapitre "gisement solaire et échanges thermiques " expose une généralité 
sur le gisement solaire et le transfert de chaleur. 
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Le deuxième chapitre est consacré à des procédés de distillation des eaux saumâtres 
par effet-serre. 

Le dernier chapitre décrire les performances thermiques et optiques du distillateur par 
des bilans énergétiques. Il décrit aussi le mode opératoire, et l’exploitation des résultats 
expérimentaux et numériques. 
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Chapitre I 


Gisement solaire et échanges thermiques 


1.1.Introduction 

La terre est considérée comme une sphère ayant un diamètre égal à 1,27 10 7 m donc 
100 lois plus petit que celui du soleil. Elle est entourée de l’atmosphère où la pression 
de 1 air diminue quant l’altitude augmente. La puissance solaire est de l’ordre de 
1 kW/m au niveau du sol. 

En ce qui concerne l’Algérie, nous avons au nord du pays 1700 kWh/m 2 /an pour 2600 
heures d’ensoleillement ; dans les hauts plateaux 1900 kWh/m 2 /an pour 3000 heures 
d ensoleillement et au sud 2200 kWh/m 2 /an pour 3400 heures d’ensoleillement. 

La dispersion fait que 1 énergie solaire est disponible, c’est une énergie qui couvre les 
besoins décentralisés et à la portée des populations rurale [12], Il est plus important pour 
notre étude de s intéresser au mouvement apparent du soleil dans la sphère céleste, les 
modes de transfert de chaleur et les rayonnements solaires et leurs appareils de mesure. 

1.2- Constitution physique du soleil[ 1,3] 

Le soleil est une sphère gazeuse composée d’Hydrogène (74 % de la masse ou 92,1 
% du volume) et d'Hélium (24 % de la masse ou 7,8 % du volume). L’hydrogène se 
transforme en hélium par fusion nucléaire, chaque seconde 700 milliards (7 1011) de 
kilogrammes d’hydrogène se transforme en 695.7 milliards de kilogrammes d’hélium. 
La diflérence -4.3 milliards de kilogrammes se transforment en énergie -38.7 10 25 
joules, c’est équivalent à 26 mille milliards de fois la consommation énergétique 
mondiale dans le même temps. Avec un diamètre de 1391000 km, le Soleil se trouve à 
une distance moyenne de 149 000 km par rapport à la Terre. A cause de l’ellipticité de 
l’orbite terrestre, la distance minimale se situe vers le 15 janvier, et la maximal vers la 
fin juin. Le Soleil n'est pas une sphère homogène, on peut y distinguer trois régions 
principales. 

a- L’intérieur: où se crée l'énergie par réaction thermonucléaire et qui est inaccessible 
aux- investigations, car le rayonnement émis dans cette région est totalement absorbé 
par les couches extérieures. La température atteint plusieurs millions et la pression un 
milliard d'atmosphère. 

b- La photosphère: l’épaisseur est d’environ 300 km qui est responsable de la presque 
totalité du rayonnement solaire que nous recevons. 
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c Chromosphère et couronne solaire: où la matière est très diluée, fait qui explique 
que bien que la température y soit très élevée (1 million de degrés) le rayonnement émis 
est très faible. 

Le soleil est la seule étoile du système solaire et la plus proche de la terre, sa lumière 
met environ 8 mn à nous atteindre. La deuxième étoile la plus proche est Proxima de 
Centaure située à 4.23 années lumières du soleil mais la grande distance qui nous sépare 
de cette étoile lait que le soleil soit la seule étoile qui assure la vie sur terre. 

Sur le plan humain, le soleil a une importance primordiale car il est à l’origine de la vie 
sur terre, en lui fournissant d’énormes quantités d’énergie, qui permet la présence de 
l’eau à l’état liquide et la photosynthèse des végétaux. Le rayonnement solaire est aussi 
responsable du climat et des phénomènes météorologiques. 

1.3- la sphère céleste[l] 

La sphère céleste est une sphère imaginaire d’un diamètre immense, avec la terre au 
centre. On considère que tous les objets visibles dans le ciel se trouvent sur la surface de 
la sphère céleste. On peut résumer les différentes caractéristiques sur la sphère elle 
même comme c’est indiqué sur (la figure 1.1). 



Figure 1.1: La sphère céleste[l] 
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1.4- Les coordonnées célestes[1,2,22] 
1.4.1 - Les coordonnées géographiques: 


Ce sont les coordonnées angulaires qui permettent le repérage d'un point sur la Terre 
(figure 1.2). 



Longitude 



90 

West 


Prime 

mprldian 


30 

East 


Figure 1.2: Les coordonnées géographiques[l] 


- La latitude <P, permet de repérer la distance angulaire d'un point quelconque à 
l'équateur. Elle est comptée positivement de 0 à 90° vers le nord et négativement de 0 à- 
90° vers le sud. 

- La longitude L, est comptée positivement vers l'est et négativement vers l'ouest, à 
partir du méridien Greenwich . 

1.4.2 - Les coordonnées horaires: 

- La déclinaison ô, c’est l’angle fait par le plan de l’équateur avec celui de 
l’écliptique. Elle varie au cours de l’année de +23°27’ à-23°27’ et détermine l’inégalité 
des durées des jours. Elle s’obtient à partir de la formule suivante: 

S = 23.45 sin x 360] Degré (1.1) 

365 

ou : 

N : le numéro du jour dans l’année compté à partir du 1 er Janvier. 
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- L’angle horaire co : il est déterminé par la rotation régulière de la terre autour de son 
axe. Il est compté positivement dans le sens rétrograde de 0° à 360° (Figure 1.3). 


Figure 1.3: Les coordonnés horaire$[22] 



1.4.3 - Les coordonnées horizontales: 

Le repérage du soleil se fait par l’intermédiaire de deux angles : 

- L’azimut a, ou l’angle que fait la projection de la direction du soleil sur le plan 
horizontal avec la direction sud. Il est compté de 0° à 360° d’ouest en est, ou de 0° à 
180° à partir du sud vers l’ouest. 

- La hauteur h, ou l’angle que fait la direction du soleil avec sa projection. Il est 
compté de 0° à 90° vers le zénith et de 0° à —90° vers le nadir. 

On appelle quelquefois distance zénithale le complément de l’angle h : z + h — 90 
(figure 1.4). 



Figure 1.4 - Les coordonnées horizontales[1,22] 


6 




Chapitre I 


Gisement solaire et échanges thermiques 



Figure 1-5: Hauteur du soleü sur site d'Adrar 


L’équation de passage des coordonnées horizontales en coordonnées horaires est la 
suivante : 

sin (a) x cos (h) = cos (Ô) x sin (to) (1.2) 

cos (a) x cos (h) = cos (S) x cos (©) x sin (O) - cos (O) x sin (ô) (1.3) 

sin (h) = cos (O) x cos (Ô) x cos (co) 4* sin (O) x sin (ô) (1.4) 

1.5- Mouvement de la terre[4,12] 

La trajectoire de la terre autour du Soleil est une ellipse dont le soleil est l’un des 
foyers. Le plan de cette ellipse est appelé l’écliptique, l’excentricité de cette ellipse est 
faible ce qui fait que 1a distance Terre-Soleil ne varie que de ±1,7% par rapport à la 
distance moyenne. 

La révolution complète de la terre au tour du soleil s effectue en une année sidérale de 
365 jours, 6 heures, 9 minutes et 10 secondes. Cette révolution faite varier les durées 
relatives du jour et de la nuit au rythme des saisons limitées par les deux équinoxes et 
les deux solstices. La terre tourne également sur elle-même autour d un axe appelé 1 axe 
des pôles en une journée sidérale de 23 heures, 56 minutes et 4 seconds. Le plan 
perpendiculaire à l’axe des pôles et passant par le centre de la terre est appelé 
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1 équateur. L axe des pôles n est pas perpendiculaire à l’écliptique, l’équateur et 
1 écliptique tont entre eux un angle appelé inclinaison et qui vaut 23°27\ 



Figurel.6- Mouvement de la terre autour de soleiI[12] 


1.6- Les temps solaires[l] 

1.6.1 - Le temps solaire vrai (T.S.V): 

Le temps solaire vrai, en un instant et un lieu donné, est l’angle horaire du soleil co. 
Il est donné sous la forme suivante : 


TSV+12+co/15 (1.5) 

co en degré. 

Si o = 0, TSV = 12h.00; Si co < 0, la matinée; Si co > 0, l’après-midi. 

1.6.2 - Le temps solaire moyen (T.S.M): 

Il est donné par la relation: 

ET = TSV - TSM (1.6) 

Avec: 

Et, représente l’équation du temps, qui corrige le T.S.V par rapport au T.S.M. 

Et= sin(2N') - 7.53 cos(N') [min] (1.7) 

N ~(N-81 ) X360/365 (1.8) 


Où : N est le numéro du jour dans l’année. 
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1.6.3 - Le temps universel (T.U): 

C’est le temps moyen du méridien de Greenwich (méridien origine). 

TU=TSM-L/15 (1.9) 

Où: 

L : est la longitude. L > 0, pour les villes situées à l’Est de Greenwich; L < 0, pour les 
villes situées à l’Ouest de Greenwich. 

TU - 'I SM, pour le méridien de Greenwich. 

1.6.4 - Le temps légal (TL): 

C’est le temps officiel d’un état, il est donné par : 

TL - TU + AH (1.10) 

Où: 

AH: le décalage horaire entre le méridien de Greenwich et l’état considéré. 

AH = 1 heure pour l’Algérie. 

L7- La constante solaire[2] 

La constante solaire est définie comme étant l'énergie totale envoyée par le soleil à 
la limite de l'atmosphère terrestre sur une surface de lm 2 placée perpendiculairement à 
la direction du rayonnement [1]. Elle est indépendante des conditions météorologiques, 
mais dépend de la distance terre-soleil. 

Cte solam =I 0 =cxT 4 x(^) 2 (1,11) 

Avec: 

R s : le rayon du soleil [mètre]; D sl : distance terre- soleil[mètre]; a : constante de 

Stephane-Boltzman: 2 7r*.k 4 /14.A 3 .<? =5,67.10 8 w/m k . 

La distance terre-soleil variant au cours de 1 année, la valeur moyenne de la 
constante solaire (I 0 ) est de l'ordre 1353W/m 2 . Cette donnée est généralement corrigée 
par la relation suivante: 

Ic=IoX(l+0.033 cos(360/365)xN ) « I 0 (l+0.034xcos(m-l)+n) (1.12) 
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Ou encore en fonction de la déclinaison: 

Ic(Ô)= I 0 x(l-(sin(6)/l 1.8)) (1.13) 

Où: ô représente la déclinaison (°), n: numéro du jour dans le mois, m : numéro du mois 
dans l'année. 

1.8- Direction des rayons solaire par rapport a un plan[2] 

Soit un plan quelconque sur la terre à un instant quelconque, repérons sa position par 
les deux angles crGt (3 tels que: 

(3 : Angle entre le plan horizontal et le plan considéré (inclinaison): 

/? >0 si le plan en face de l'équateur, (3 =0 pour le plan horizontal , f3 = 90° pour un 
plan vertical. 

a : Angle entre le méridien local et la normale au plan 

a ~0 vers le sud, a <0 vers l'est, a>0 vers l'ouest, a-n vers le nord. 

Soit "i" l'angle entre le rayon incident et la normale au plan considéré. 

Pour une surface orientée vers le sud et inclinée d'un angle (3 sur l'horizontale, La valeur 
de "i" est donnée par la relation suivante : 

Cos(i) =Cos(œ) xCos(ô ) x Cos(d> -/?) +Sin{ô) xSin(<P - /3) (1.14) 

1.9- Eclairement hors terrestre[4,21] 

Le rayonnement hors atmosphère couvre un très large domaine de longueur d ondes il 
ne dépend d'aucun paramètre météorologique, mais il est en fonction de quelque 
paramètre astronomique et géographique tel que: la latitude du Heu (p, la déclinaison 
solaire 5 et l'angle horaire au coucher du soleil o c . Sur une surface horizontale, et 
pour le jour n, le rayonnement extraterrestre G D ( MJ/m .jour) est obtenu à laide de 
l'équation suivante: 


G o - 


24 x 60 


K 



r ( 2xn\\ 

coc.Sin (<p ).Sin (S ) + 

.Gsc . 

1 + 0 034 M 365 Jj 

Cos ( (p).Cos (S).Sin (ûjc) 


(1.15) 
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• L’Indice de clarté 

C est le rapport entre le rayonnement horizontal au sol et le rayonnement extraterrestre, 
notée Kt est défini par l’expression : 

K T = -p- (1.16) 

Go 

Go :c’est le rayonnement horizontal au sol • 

• L’éclairement 

L éclairement est la puissance de rayonnement reçu par unité de surface. 
L’atmosphère ne transmet pas au sol la totalité du rayonnement solaire. Les nuages, les 
gouttes d eau, les molécules de gaz, les particules solides en suspension dans 
1 atmosphère empêchent la totalité du rayonnement solaire à atteindre le sol. 

1.10- La radiation solaire au niveau du sol [2] 

Le rayonnement solaire subit un certain nombre d’altérations et d'atténuations 
aléatoires en traversant l'atmosphère; réflexion sur les couches atmosphériques, 
absorption moléculaire, diffusion moléculaire et par les aérosols (la poussière, les 
gouttelettes....). Au niveau du sol, à cause de la diffusion, une partie du rayonnement 
est diffuse (c’est-à-dire isotrope). L'autre partie dite directe est anisotrope. Le 
rayonnement global est par conséquent la somme du rayonnement direct et du 
rayonnement diffus. 

1.10.1- Composantes du rayonnement solaire: 

L'atmosphère ne transmet pas au sol la totalité uu rayonnement solaire qu’elle reçoit. 
Les composantes du rayonnement solaire incident instantané sont: 
a- Le rayonnement direct, provenant du disque solaire et qui subit l’atténuation due à 
la traversée de l'atmosphère. L'expression du rayonnement direct pour un plan orienté 
vers le sud est donnée par l'expression: 

/ D = / c x a x exp(-^~x P/1000).cosi (1.17) 

Avec: a et b: Coefficients traduisant les troubles atmosphériques. 
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I ,c tableau (1.1 ) rassemble quelques valeurs de a et b. 

Sinh = Cos(5) x Cos(üj) x Cos(0) + Sin(S ) x Sin(0) (1.18) 
h: Etant la hauteur du soleil. 

ù), Ô, 0 Sont respectivement l'angle horaire, la déclinaison, et la latitude [°]. 

La pression atmosphérique "p" en fonction de l’altitude est illustrée par le tableau (1.2) 


Tableau 1.1: Valeurs des coefficients de trouble atmosphérique 



Ciel pur 

Condition normale 

Zones industrielles 

A 

0.87 

0.88 

0.91 

B 

0.17 

0.26 

0.43 


Tableau 1.2: Variation de la pression atmosphérique avec l'attitude 


Attitude 

0 

500 

1000 

1500 

2000 

5000 

P (mbar) 

1000 

950 

900 

850 

800 

500 


b- Le rayonnement diffus: provient de touts les éléments irradiés par le soleil. 
Le flux diffus reçu par la surface du capteur se caractérise par deux composantes du 
flux à savoir: 

Le flux diffus émit par la voûte céleste (q>i), donné par la relation: 

cp i=(p 0 x ( 1 +cos(f3)/2) (1.19) 

Le flux émis par le sol ( 92 ) : 

(p2=px(l-cos(p)/2)xI c XT dir xsin(h)+ (p 0 (1.20) 

tpo= Icxsin(h) x(0.271-0.2939xx dir ) (1.21) 

avec: 

Tdir : La transmissivité totale de la couche atmosphérique dont l'expression est : 

x air = a x exp(—^xp/100) (1.22) 

p L'albédo du sol donnée par le tableau (1.3) 


Tableau 1.3: Albédo de différents sols 
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Nature du sol 

L’albédo 

Herbe fraiche 

0.20 

Herbe sèche 

0.30 

Terra cultivée, nue 

0.16 

Gravier 

0.22 

Neige fraiche 

0.85 

Etndu d’eau 

0.05 


Le flux diffus total s'écrit par conséquent de la manière suivante: 

Id~ <Pi+ <P2 (L23) 

c- le rayonnement global: Le rayonnement global sur une surface quelconque au 
Voisinage du sol est la somme de ces deux composantes: 

G = I D + I d (L24) 


1.11- Propriétés optiques du vitrage 

1.11.1- Facteur solaire: 

Le facteur solaire (FS) d’un vitrage est le rapport entre l’énergie totale entrant dans 
le local à travers ce vitrage et l’énergie solaire incidente. Cette énergie totale est la 
somme de l’énergie solaire entrant par transmission directe, et de l’énergie cédée par le 
vitrage à l’intérieur à la suite de son échauffement par absorption (flux transmis + flux 
absorbé réémis vers l’intérieur)(figurel.7). 

incident Flux réémis vers 



Figure 1.7: Propriétés optiques du verre[4] 
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1.11.2- Effet de serre: 

Lorsque le rayonnement solaire frappe une paroi vitrée, une partie sera réfléchie, une 
partie sera absorbée puis réémise et le reste sera directement transmis. Ce rayonnement 
direct transmis à l’intérieur des locaux est partiellement absorbé par les parois qui 
s’échauffent et par la suite le rayonnement infrarouge est réémis dans toutes les 
directions. Le vitrage étant pratiquement opaque au rayonnement réémis, celui-ci reste 
piégé à l’intérieur du local, y entraînant une augmentation de la température. C’est ce 
qu’on appelle l’effet de serre figure(L8). 



L12-Echanges thermique[3] 

Le transfert de chaleur est l'un des modes d'échange d énergie les plus connus. Le flux 
de chaleur entre deux corps est dirigé spontanément du corps froid. 

Un corps quelconque chauffé absorbant une quantité de chaleur (Q) peut. 

• Soit voir sa température augmenter d'une valeur (A7),ou 

• Soit conserver sa température Tq — Cte en changeant détat (du solide vers le 
liquide, du liquide vers la vapeur, du solide vers la vapeur). On dit alors qu’il 
échange une quantité chaleur dite chaleur latente de changement d état. 

1.12.1- Echange de chaleur par conduction: 

La conduction est définie comme étant le mode de transfert de chaleur provoqué par 
la différence de température entre deux régions d’un milieu solide, liquide ou gazeux. 
L'effet macroscopique observable est une égalisation des températures du système. 
Cependant si certaines zones sont maintenues à température constante par apport de 
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chaleur (réservoir de chaleur) ou évacuation de chaleur (puits de chaleur), il s'établit un 
transfert continu de la chaleur de la région chaude vers la région froide. 

d Q = S ' dt (>- 25 ) 

dQ: énergie élémentaire en Joule; X: Conductivité thermique du matériau en W.m^.K 1 ; 
S: section en m 2 ; dt: temps élémentaire en s; gradient de température en x [K.m 1 ]. 

1.12.2- Echange de chaleur par convection: 

Dans un processus de convection le transfert de chaleur entre un solide et un fluide 
l'énergie étant transmis par déplacement du fluide. Ce mécanisme de transfert est régi 
par la loi de Newton: 

Q = h.S(T p -T„) (1.26) 

Q: Flux de chaleur transmis par convection [W]; 
h : coefficient de transfert de chaleur par convection [W/m 2 .°C]; 

T p : Température de la surface solide[°C]; 

Too : Température du fluide de la surface solide [°C]; 

5:Aire la surface contact solide/fluide [m ]. 

Selon le mécanisme qui génère le mouvement du fluide, on distingue . 

• La convection naturelle ou libre 

• La convection forcée 


1.12.3- Echange de chaleur par rayonnement: 

Tout corps chauffé émet par sa surface extérieure un rayonnement électromagnétique 

dont la puissance est fonction de sa température ; C'est le rayonnement thermique. 
L'analyse spectrale de ce rayonnement montre une majorité de courtes longueurs d'onde 
aux très hautes températures et une majorité de grandes longueurs d'onde pour des 
températures inférieures à 500 K. 

D'après la loi de Stéphan, le flux de chaleur échangé entre la surface S et le milieu 

ambiant peut s'écrire : 

L . /T 
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Glapis ~ 1 


o : constante de Stéphan Boltzmann 5.67 10' 8 [W.m^.K" 4 ]; 

5 : surface d'échange |m 2 ] ; 

6 : coefficient d'émission de la surface (s=l pour un corps noir, e«l corps brillant); 
T s : température de surface du corps [k]; 

T a : température ambiante[k]. 


1,12.4- Echange de chaleur par changement de phase: 

Lors de changement de phase d'un corps un échange de chaleur s’effectue, le flux 
chaleur échangé est: . 

Q-m.L (1.28) 


Tel que: 

m: la masse de corps 

L: La chaleur latente massique de changement d'état représente la quantité de chaleur 
nécessaire pour permettre le changement d'état d'une unité de masse (ou une mole) d un 
corps pur d'une substance sous une pression constante. Elles s expriment en J.kg 
(J.mol 1 ) dans le système légal. Plus couramment on trouve les unités kj.kg (kJ.mol ). 

Suivant la convention thermodynamique, les chaleurs latentes de changement détat 
sont positives si on doit fournir de l’énergie à un corps ou négatives si le corps libère de 
l’énergie. 
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1.13-1 jCs appareils de mesure des rayonnements solaire 
1.13.1- Pyranomètrel 14J: 

Le pyranomètre mesure le rayonnement sur une surface plane (Watt/m ). Le 
rayonnement mesuré représente la somme du rayonnement direct produit par le soleil et 
du rayonnement diffus. Ils ne nécessitent aucune alimentation extérieure. Ils génèrent 
une tension qui, typiquement, équivaut à : 10(mV/(kW. m 2 ). Chaque pyranomètre est 
calibré individuellement conformément à WRR (World Radiometric Référence). 



Figure I.9.a : Un'pyranomètre pou7ÎT~ Figure I.9.b : Un pyranomètre bande pare-soleil 
^ i 1 » musnro ri» ravnnnement sohirc dinuSi 


mesure du rayonnement solaire 
global[141 


pour la mesure du rayonnement solaire diffus, 
où l’ombre de la bande apparait sur la demi 
sphère de P instrument [14]_ 


1.13.2- Le pyrhéliomètre[15|: 

Le pyrhéliomètre est conçu pour effectuer des mesures du rayonnement solaire direct 
en incidence normale, sans surveillance. Un pyrhéliomètre est un instrument conçu 
spécialement pour mesurer le flux lumineux solaire du faisceau direct avec un champ 
visuel limité à 5 ». Ce résultat est obtenu grâce à la configuration du tube collimateur, 
doté d’ouvertures de précision, et à la conception du détecteur. L’ouverture avant est 
équipée d’une vitre de quartz pour protéger l’instrument et pour faire office de filtre 
faisant passer le rayonnement solaire à des longueurs d’ondes comprises entre 200 nm 
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1.13.3- Hcli<>Rraphe[ 16] : 

Il est utilise pour mesurer le nombre d'heures d'ensoleillement dans une journée. Une 
boule de cristal maintenue par un support amovible (qui peut bouger) constitue la base 
de l'héliographe. Derrière la boule de cristal est placé un carton sur lequel les mesures 
s'inscriront. L'appareil est fixé sur un support d'une hauteur approximative de 1,5 m. 
Chaque soir, après le coucher du soleil, il faut enlever le carton qui a été brûlé par les 
rayons du soleil ayant traversé la boule de cristal. Celle-ci agit comme une loupe en 
concentrant les rayons sur le carton. Le déplacement du soleil fait changer le point de 
convergence des rayons sur le carton, ce qui produit une ligne brûlée sur celui-ci. Si un 
nuage cache les rayons du soleil, il y aura sur le carton un endroit, correspondant à cette 
heure, qui ne sera pas brûlé. À l'aide d'une règle graduée en heures, on mesure la 
longueur des brûlures sur le carton. 

Ainsi, en additionnant le temps de toutes les sections, il est possible d obtenir le nombre 
d'heures d’ensoleillement total de la journée en tenant compte de la présence ou de 
l'absence de nuages. Cet appareil fut inventé par J.F. Campbell vers 1853 et modifié 
vers 1879 par G.G. Stokes. On doit placer cet appareil loin des édifices et des choses qui 
pourraient produire une ombre sur lui. 



Figurel.il : Héliographe[16] 


1.14- Conclusion: 

Ce chapitre a permet de décrire les paramètres astronomiques et radiatifs . l'effet de 
serre, les composantes de rayonnement solaire et leur instiuments de mesure utilises. 
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II. 1- Introduction 

Le dessalement est le processus de traitement d’eau qui consiste à enlever les sels 
de l'eau saline. Cette technologie a été réellement utilisée pour la première fois il y a 
environ soixante ans, durant la deuxième guerre mondiale, quand les bases militaires 
dans les régions arides ont eu besoin d’eau pour leurs troupes. Le gouvernement 
américain a investi beaucoup pour le développement des technologies du dessalement. 


Par la suite, des stations de dessalement de capacité 8000 (m 3 /j) ont été installées dans 
diverses régions du monde. 

Elles sont utilisées spécialement pour le dessalement de l’eau de mer. En 1970, de 
nouvelles technologies, des processus commerciaux membranaires comme 
l’électrodialyse (ED) et l’Osmose inverse (RO), sont apparues et commençaient à être 
utilisées. 1/HI) est appropriée pour dessaler l’eau saumâtre parce qu’elle est beaucoup 
plus économique que la distillation alors que le RO est utilisée pour dessaler à la fois 
l’eau saumâtre et les eaux de mer qui sont très salées. 

Dix ans plus tard, le dessalement est devenu une entreprise commerciale mature. Un 
grand nombre de chercheurs ont travaillé dans ce domaine pour améliorer le rendement 
de ces technologies et atteindre une efficacité meilleure. A partir de 1990, 1 utilisation 
des technologies de dessalement est devenue une chose banale. Ces technologies étaient 
utilisées pour produire de l’eau potable pour préserver la vie dans certaines régions 
comme les pays de la région aride et les pays de Moyen-Orient. Par la suite, des 
technologies diverses ont été développées et certaines d’entre elles sont devenues de 

plus en plus efficacités^12]. 


II.2- Les différents procédés de dessalement[3,12,13,22] 

Le dessalement de l’eau de mer et des eaux saumâtres, par les différents procédés 
existant, connaît un développement important. Parmi les plus connus, il y a les 


processus de distillation thermiques comme . 

• La distillation à simple efiet (simple effet distillation [SED]) 

• La distillation à multiples effets (multiple effet distillation [ME ]) 

• La distillation à multiples étages (multiple stage flash [MSF]) 

• La distillation par compression de la vapeur (vapeur compression [V ]) 
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Les procédés Membranes sont les procédés dont la source principale d’énergie est 
réieetricitc comme : L osmose inverse (reverse osmosis [RO]) et L’électrodialyse [ED]. 
De plus il > ad autres procédés comme le procédé d’échange d’ions, et le procédé de 
congélation. 

a- La distillation à simples effets [SED) 

Une circulation de vapeur permet de chauffer dans un évaporateur l’eau salée, qui 
s’évapore sous une pression légèrement supérieure à la pression atmosphérique, on peut 
alors condenser directement la vapeur sortant de l’évaporateur, sans récupération de la 
chaleur de condensation, ou bien, au contraire, on peut envoyer la vapeur dans un 
condenseur en vue du préchauffage de l’eau salée qui arrive dans l’installation. 



Figure II.1- Distillateur simple effet[22] 


b- La distillation à multiples effets [MED] 

Pour augmenter le rendement de l’opération, il y a intérêt à récupérer la chaleur 
de condensation et à l’utiliser comme chaleur d apport pour une nouvelle distillation, on 
réalise alors plusieurs effets en cascade, la chaleur résiduelle est récupérée et utilisée 
pour le chauffage de l’eau brute à traiter. Dans les industries de traitement divers, on 
réunit ainsi jusqu'à six ou sept évaporateurs - condenseurs (le premier évaporateur et le 
dernier condenseur ayant une seule fonction), et peut être combinée avec une centrale 
électrique thermique. Du point de vue technologique deux types peuvent etre utilisée a 
long tube vertical (LTV) inventé en 1905 par le français KESTNER ou à tube 
horizontal. Pour fixer le nombre d’effet, on tient compte de la différence globale de 
température disponible entre la source chaude et froide, et du gradient de température 
nécessaire dans chaque évaporateur pour assurer efficacement la transmission de 
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chaleur de la vapeur de chauffe à l’eau saline à distillée. Il dépend aussi du prix de 
revient de la vapeur d’alimentation, donc du coût de combustible employé. 



Figure 11.2: Distillation à multi effets (M.E.D)[22| 

c- La distillation avec vaporisation instantanée à étages multiples ou par 
détentes successives [MSF] 

L’installation comprend un certain nombre d’enceintes ou cellules, thermiquement 
isolées disposées en série dans lesquelles se produit la vaporisation instantanée par 
détente. 

L’eau saline à traiter est chauffée sous pression, sans évaporation, puis introduit dans les 
cellules sous vide partiel ou elle circule successivement, a chaque étage la vapeur 
produit par détente et va se condenser a la partie supérieur au contact d’un faisceau 
tubulaire dans lequel l’eau d’alimentation circule. 



Figure II.3- Distillateur par détentes successives[22] 
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d- La distillation avec compression de la vapeur (thermo compression)[VC] 
1/eau salée arrive dans un évaporateur, ou elle se vaporise aux environs de 100°C grâce 
à un serpentin auxiliaire de vapeur qui est ensuite mis hors circuit. La vapeur d’eau à la 
pression atmosphérique et vers 100°C est alors aspirée par un compresseur qui la porte à 
une pression relative de l’ordre de 0,2 bar, avec une surchauffe de l’ordre de seulement 
6°C, la vapeur surchauffée passe ensuite dans le serpentin principal de l’évaporateur où 
elle se condense. Les avantages de la thermo compression sont son haut rendement, on 
peut retirer environ 0,8 m 3 d’eau douce de 1 m 3 d’eau brute. 


Compresseur 
de vapeur 



Figure IL4 : La distillation avec compression de la vapeur [131 


e- L’osmose inverse [RO] 

L'osmose est le transfert de solvant à travers une membrane sous l'effet d’un 
gradient de concentration. Si on considère un système à deux compartiments séparés par 
une membrane semi sélective et contenant deux solutions de concentrations différentes, 
l'osmose se traduit par un flux d'eau dirigé de la solution diluée vers la solution 

concentrée. 


Si on applique une pression sur la solution concentrée, la quantité d'eau transférée par 
asmose va diminuer. Avec une pression suffisamment forte, le flux d'eau va meme 

s’annuler: cette pression est nommée la pression osmotique (en faisant l'hypothese que 

, x o- ccp la valeur de la pression osmotique, on 

la solution diluée est de l'eau pure). Si on dépas 

Hn flux osmotique: c'est le phénomène 

observe un flux d'eau dirigé en sens inverse 


d'osmose inverse. 
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OSMOSE OSMOSE INVERSE 


Figure II.5: Les techniques de dessalement par osmose[3] 
f- L’électrodialyse [ED] 

Un élcctro-dialyseur est formé essentiellement d’une cuve dans laquelle on place un 
certain nombre de membranes sélectives, en faisant alterner les membranes anionique et 
les membranes cationiques, qui délimitent une suite de compartiments .On fait circuler 
dans cette cuve l'eau saline, tandis qu’on la fait traverser par un courant électrique 
continu grâce à des électrodes, donc on va voir une migration des ions de l’eau saline à 
travers ces membranes imperméables à l’eau. 

Une membrane sélective est infranchissable par les ions d’une catégorie tels que les 
cations, mais laisse pénétrer sur une face de la membrane les ions de l’autre catégorie 
(anions) et ne laisse pas pénétrer sur l’autre face, donc la membrane était perméable au 
anions et imperméable aux cations, on dit que c‘est une membrane anionique et on 
définirait de façon analogique une membrane cationique, possédant la propriété inverse. 



Eîill douce 

Figure [1.6: Schéma (le principe de rélectrodialyse[13] 
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g- Echange d’ions 

Le phénomène d'échange d'ions se produit quand on met en contact une eau 
saumâtre avec certaines substances insolubles ayant la propriété d’échanger certains de 
leurs ions avec les sels de 1 eau, les produits employés actuellement sont des produits 
insolubles dénommés échangeurs d'ions (résines synthétiques sous la forme de graines 
d environ 0,3 à 1,5 mm). Cette méthode est applicable seulement pour les eaux 
contenant jusqu'à 8 et même 10 g/L. 


h- La méthode de congélation 

La température de congélation est presque indépendante de la pression, elle est 
déterminée par la concentration saline et varie depuis 0°C par l’eau pure jusqu'à environ 
-2°C pour l’eau de mer et -4°C pour une concentration double de celle de l’eau de mer, 
le dessalement par congélation est basé sur le fait qu’une fraction d’une masse donnée 
d’eau saline se dérobe à la congélation en s’enrichissant en sels, la partie congelée étant 
formée d’eau pure. 

Quelques remarques : 

1) D’une façon générale les procèdes thermiques représentent, environ 75% du volume 
d’eau dessalée, l’osmose inverse occupe environ 22% du marché mondial le reste est 
porté sur le compte des autres procédés. 

2) L’utilisation de toutes ces procédées exige, cependant, de gros investissements que la 
plupart des pays en voie de développement ne possèdent pas, ces régions ont besoin de 
systèmes de dessalement performant certes, mais leur capacité de production doit être 
appropriée à la consommation modeste des habitants et aux ressources énergétiques 
naturelles qu’elles possèdent sur leur territoire, par exemple, 1 énergie solaire[12]. 
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11.3- Distillation solaireI4,13,18| 

La distillation solaire est une technique qui utilise le rayonnement solaire pour 
chauffer 1 eau salée ou saumâtre dans un bac couvert par une vitre inclinée. Une partie 
de l'eau s’évapore et la vapeur produite se liquéfie sur la surface intérieure de la vitre. 
La distillation solaire a été utilisée depuis plusieurs siècles; mais la plus grande 
réalisation est apparue en 1872 prés de Las-Slinas dans le nord du chili. Le modèle a été 
fabriqué par Carlos Wilson, un ingénieur Suédois avec une surface vitrée de 5000m 2 , 
produisant jusqu'à 20 m 3 d’eau douce par jour, à partir d'une eau saline contenant 140g/l. 
La distillation solaire connut une période creuse entre 1880 et la première guerre 
mondiale. Kn 1920, Kaush utilisa des réflecteurs en métal pour concentrer les rayons 
solaires. Une installation due à Ginistous groupant 30 appareils de 1m 2 fut réalisée et 
exploitée en Tunisie en 1929. Au début des aimées 30, Tifïnov proposa un distillateur 
incliné. Maria Telles inventa en 1945, un nouveau type de distillateur solaire de forme 
sphérique. Kn 1969, Cooper proposa une simulation pour analyser les performances d'un 
distillateur solaire type serre. Depuis les années 70, plusieurs autres types de 
distillateurs ont été élaborés et étudiés, parmi lesquels: 

* Le distillateur à effet multiple 

■ Le plateau incliné ou distillateur solaire incliné 

■ Le distillateur solaire type mèche ou type multiple mèche. 

11.4- Principe de la distillation solaire[2] 

Les distillateurs solaires utilisent tous le même principe de fonctionnement, 
néanmoins, ils peuvent différer les uns des autres en conception et en matériau. Le 
fonctionnement d'un distillateur solaire repose sur l'effet de serre. 

L'énergie solaire entrant par la couverture est absorbée par l'eau et le fond noirci, puis 
elle est transmise à l'air par conduction, convection et rayonnement. Comme les parois 
sont isolées et le verre opaque aux radiations infrarouges de grandes longueurs d'ondes 
l'énergie solaire se trouve ainsi emprisonnée (effet serre) et 1 air chauffé (figure 1 ) 

La température de la couverture en contacte avec 1 air atmosphérique étant infé 
celle du mélange air vapeur d'eau, il se produit une condensation de la vapeur d’eau sur 

ta surface intérieure de la vitre. 




( 11 1 Mc disüllatiü,, de., c.» ux ^mâlr» „ a r effet-.,,-.... 



Figure II.7: Principe de fonctionnement d'un distillateur solaire à effet de serre[2J 


11*5- Les différents types de distillateurs solaires [7,13,18] 

Les distillateurs solaires connus se répartissent en deux grandes groupes, ceux dont 
la surface de condensation est également la couverture transparente au rayonnement 
solaire (mode direct), et ceux dont la surface de condensation est séparée du chauffage 
de l’eau (mode indirect). 

11.5.1- Distillation solaire à effet de serre: 

Soit une serre fermée et exposée au soleil, à T intérieur de laquelle se trouve une lame 
d’eau de mer ou d’eau saumâtre de quelques centimètres d’épaisseur. L’air à 1 intérieur 
de la serre est surchauffé et saturé de vapeur d’eau (douce) qui se condense au contact 
de la paroi en verre relativement froide. Les gouttes d’eau douce peuvent être recueilles 
en bas du vitrage dans une gouttière (figure II.8). 



Figure IL8: Distillation solaire à effet de serre[13] 
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Il 5.2- Distillateur solaire incliné à cascades: 

Cet appareil présente deux avantages, la lame d'eau est très faible et l’orientation par 
rapport au rayonnement incident se rapproche davantage de l’optimum. Un modèle de 
ce type (figure 11.9) n’a été exécuté qu’en petite taille. Ces appareils fonctionnent bien 
mais leur construction et leur entretien sont coûteux. 



Figure 11.9-Distillateur incliné à cascades (ou à marches)[7] 


11.53- Distillateur solaire incliné à mèche: 

Il se compose d’une très mince couche de textile à mailles fines, imbibée d’eau qui se 
maintient d’elle même au contact de la plaque jouant le rôle d évaporateur. Les forces 
de tension interfaciales, nettement supérieures à celles de la pesanteur, jouent un rôle 
très important dans la formation du film capillaire, de faible inertie thermique, en 
empêchant toute formation de bulles d’air. Ce distillateur est composé d'un capteur - 
évaporateur, d'un condenseur et d'une alimentation. La première cellule p 

rayonnement solaire qui passe à travers la couverture. La vapeur se cond 
paroi opposée et la chaleur dégagée par cette condensation permet l’évaporation du film 

qui ruisselle sur l’autre face de cette meme paroi. 

1.5.4- Distillateur solaire sphérique: 

C'est un distillateur eu plereigta «sp™»- * * , 

j Tt 1 AX T e rayonnement incident transmis par 1 

condensation demi- sphérique (ligure II. 10). L y 

« - ,p M re supérieur, p* dans „ - e, .ed, .- £ - ‘ ~ 
dans celui-ci. Une partie de l’eau s évaporé, et la vap inférieure à travers 

la face intérieure de la vitre. Le distillât passe dans le dem P 
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l'espace annulaire entre le bac et la sphère. Le balayage de la surface de condensation 
est réalisé au moyen d'un essuie - glace entraîné par un moteur. L’avantage de ce 
système de balayage est de maintenir la surface de condensation, constamment 
transparente au rayonnement, et d'assurer un drainage rapide des gouttelettes. 



Figure IU0- Le distillateur solaire sphérique[2] 


11,5.5- Distillateur solaire à étages multiples avec récupération de chaleur: 

Seul, le bassin inférieur est alimenté en énergie thermique (figure LU). Lorsque l’eau 
de ce bassin est réchauffée, elle s’évapore et se condense à la surface inferieure du 
bassin conique situé au-dessus. 11 y a une formation de gouttes d’eau qui s’écoulent vers 
une gouttière qui les conduit au stockage. A la condensation de la vapeur d’eau, 
l’énergie de conversion de phase, appelée aussi l’enthalpie d’évaporation, est libérée. 
Cette énergie réchauffe l’eau contenue dans l’étage supérieur qui s’évapore à son tour et 
se condense à la surface inférieure de l’étage suivant et ainsi de suite...De cette maniéré, 
l’énergie, prévue pour réchauffer l'eau du bassin inférieur, est récupérée, du moms en 
grande partie, pour être plusteurs fois utilisée. Ce processus de récupération de chaleur a 
un elfe, multiplicateur sur la production d’eau potable. Le remplissage du délateur 
s’effectue à partir du dernier étage, le système est alimente en eau saline a un de ît 
fois plus élevé que celui de l’eau évaporée. Cette quantité d’eau mjeetee 

continuellement à travers les tuyaux de débordement a pou 

- de compenser la quantité d’eau évaporée, ^ 

. , | ç bassins et de contribuer a conserver 

- d’éviter une accumulation des sels dan 

les gradients de température entre les étages. ^inhnn- il 

. . . r«filiation fonctionne en thermosiphon, u 
Le circuit entre le capteur solaire à 1 unité e 1 . p r éviter 

^ HVnereie supplémentaire. Pour éviter 
n'y a donc pas de pompe nécessitant une sourc 
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une surpression, une huile végétale est utilisée comme fluide caloporteur car la 
température peut dépasser les 100°C. 



Figure II.Il: Schéma d’un distillateur solaire à étages multiples[13] 


II.6- Paramètres et caractéristiques de fonctionnement:[8] 

II.6.1- Les paramètres influants sur le fonctionnement du distillateur: 

Ont peut classer les paramètres qui influent sur le fonctionnement du distillateur 
comme suit: 

- Les paramètres externes. 

- Les paramètres internes. 


n.6.1.1- Paramètres Internes: 

* Paramètres de construction: 

- Couverture: 

La couverture verrière intervient essentiellement par sa nature qui peut etre en 
verre ou en plastique. Elle a deux rôles à jouer : d’une part c est un filtre sé 
rayonnement solaire (le rayonnement émis ou réémis par 1 intérieur du 
température basse qui ne la traverse pas) et d’autre part, c est une surfa 
condensation de la vapeur d’eau : une bonne mouillabüité est nécessaire pour éviter la 

condensation en gouttelettes qui ont tendance à retomber dans 
réfléchissent une partie importante du rayonnement incident. Pour assurer 

ruissellement de l’eau condensée vers la gouttière (collect ) 
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- L’inclinaison: 

Son inclinaison par rapport à l’horizontale, pour déterminer la quantité d’énergie 
solaire introduite dans le distillateur et pour minimiser la distance entre la saumure et la 
vitre, l’angle d’inclinaison doit faire l’objet d’un choix judicieux. 

L’inclinaison influe aussi sur les équations des bilans énergétiques des différents 
constituants du distillateur. Elle dépend du fonctionnement du distillateur durant 
l’année: 

• Fonctionnement estival: fi — (p 

■ Fonctionnement hivernal: fi — <P 

■ Fonctionnement annuel: p= <P -10° 

- L’absorbeur: 

Les études faites dans ce domaine montrent que la surface, absorbante peut être 
construite de plusieurs matériaux (bois, métal, béton, matière synthétique ou en verre 
ordinaire). Le choix de la matière de la surface absorbante ou du bac noir dépend de son 
inertie thermique, de la résistance à l’oxydation par l’eau et des dépôts minéraux. 

* La distance entre la surface d’évaporation et la sur face de condensation: 

STCUNANTHAN et HANSEN ont montré que la performance du distillateur augmente 
quand la distance entre la saumure et la vitre diminue. 

•L’hauteur du niveau de la masse d’eau à distiller: 

L’épaisseur de l'eau saumâtre joue un rôle très important. La production autant plus 
importante pour un distillateur à faible épaisseur de saumure, mais pour un distillateur 
de grande épaisseur de la saumure, le maximum de la production n est observé que peu 
après le coucher de soleil. 

'L’isolation des faces latérales du distillateur: 

L’isolation des faces latérales est pour but de éliminées les pertes (transfert de chaleur 
avec le milieu extérieur). 


v Paramètres thermophysiques: 

Les paramètres thermophysiques sont représentés par : 

La couverture et la surface absorbante : la conductivité thermique, la chaleur spé ’ q 
et la diffusivité thermique. 

** saumure: la conductivité thermique, la chaleur spécifique, la viscosité dynamique, la 
viscosité cinématique, la chaleur latente de vaporisation, le coefficient d 
dermique et la diffusivité thermique. 
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U mélange air-vapeur : la conductivité thermique, la chaleur spécifique, la viscosité 
dynamique, la viscosité cinématique, le coefficient de dilatation thermique. 


<• Paramètres optiques: 

Ce sont les paramètres optiques (l’émissivité, l’absorptivité, la réflectivité et la 
transitivité) de la sur face de la couverture et de la surface absorbante. 

$ Paramètres de fonctionnement: 

Ces paramètres sont les différentes températures du système et du fluide, en particulier 
la température de l'eau à distiller. 

MU- Paramètres Externes: 

Les études faites par plusieurs chercheurs ont montré que les paramètres les plus 
fréquemment retenus sont: 

4 Paramètres météorologiques 
i) Intensité du Rayonnement Solaire 

b) la température ambiante de l'air et la vitesse du vent 

Aux températures de saumure élevées, augmenter la différence de température entre la 
saumure et la vitre en augmentant la vitesse de vent, en baissant la température 
ambiante. Par contre, aux températures faibles de la saumure, et en augmentent la 
vitesse du vent revient à faire chuter le taux d’évaporation. 

c) L’intermittence des nuages : 

La quantité d'énergie solaire reçue par une surface dépend des paramètres A, et B qui 
! ^pendent de la nature de ciel. Cela revient à dire que le passage des nuages diminue 
1 absorptance effective moyenne. 

* Paramètres Géographiques 

Les paramètres géographiques tels que la longitude, la latitude et la hauteur d 
sont à prendre en considération dans cette étude. 


I 
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II.7- Caractéristiques du distillateur[3] 

Plusieurs grandeurs sont définies afin de caractériser la production d’eau d’un 
distillateur. On distingue constamment, le rendement, l’efficacité interne et globale, et la 
perfoimance. Le taux de production ; 

• Le rendement: 

C’est la quantité d’eau produite par unité de surface de plan noir et par jour. 
L’inconvénient majeur de ce critère est qu’il ne fait pas mention de l’énergie solaire qui 
anive sur le distillateur. 


à : Débit massique de distillât. 

4: Chaleur latente de vaporisation. 

G; Puissance de rayonnement solaire incident par unité de surface. 


' L'efficacité: 

* L'efficacité globale: 

L’efficacité globale est définie par le rapport: 



(H-2) 


L'efficacité interne: 

Si 1 efficacité globale rapporte la quantité d’eau produite à l’énergie solaire tombant sur 
me surface horizontale, elle ne fait pas mention de la quantité d’énergie entrant 
Vilement dans le distillateur, d’autant plus qu’un distillateur est construit en général 
pour un lieu donné et avec une pente de couverture fixe. Un changement de pente et de 
la localité fait varier la quantité d’énergie qui entre dans le distillateur. On définit donc 
** c fficacité dite interne qui tient 
compte de tous ces paramètres et qui est définie par: 

„ _ 2s (H-3) 

m ~ Qe 
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« La performance 

u souci de caractériser un distillateur d’une manière absolue a amené ainsi 
Satcunanathan et Hansen |9| à définir le facteur de perfotmance brut (F.P.B) et le 
Acteur de performance horaire (KP.H): 

p p g _ Quantité d'eau produite au bout de n h 

Quantité d’énergie entrée au bout de T <h (^-4) 


ppg _ Quantité d'eau produite au bout d'une heure 
Quantité d'énergieentrée au bout d'une heure 

E8- Technologie de la cou vertu re[7] 

| Dans la plupart des appareils, la couverture est réalisée soit en verre soit en plastique. 
Elles a deux rôles à jouer: d’une part c'est un filtre sélectif de rayonnement solairefle 
nyonnement solaire la traverse, le rayonnement réémis par l’intérieure du bassin à 
température basse qui ne la traverse pas} et d’autre part, c’est une surface de 
condensation de la vapeur d’eau : une bonne mouillabilité est nécessaire pour éviter la 
condensation en gouttelettes qui ont tendance à retomber dans le bassin et qui 
réfléchissent une partie importante du rayonnement incident. Pour assurer le 
ruissellement de l’eau condensée vers la gouttière (collecteur), le toit est en pente (5 à 
20°) en général). L’emploi du verre permet de satisfaire ces exigences ; le seul 
inconvénient étant son prix relativement élevé. C’est pourquoi différents auteurs 
préconisent l’emploi de films plastique (Chlorure de polyvinyle, polyéthylène, nylon, 
Iddar). L’expérience a montré que les films plastiques sont fragiles, vieillissent trop vite 
résistent mal aux ouragans ! Toutefois, l’emploi du plastique est systématique pour 
les petits distillateurs portatifs ou embarqués à bord de canots de sauvage. 


^Technologie du bassin[7J 

bassin est réalisé soit en béton dans le cas de grands distillateurs, soit en métal 
011 en polyester stratifié dans le cas des petits distillateurs. L’isolation thermique entre 


,eau ^ le terrain d’assise peut être réalisée soit par un revêtement intérieur au bassin, 
« masse » par l’emploi de ciment d’amiante. L’étancheite du fond du bassin 

/la 


dans la 
esî balisée 


par un tapis de revêtement en 


,-nr.nichouc synthétique ou en toile de verre 
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bitumée. Ce même revêtement de couleur sombre absorbe le rayonnement solaire. 
Eventuellement F utilisation de l'énergie absorbée peut être améliorée par la mise en 
place dans l'eau du bassin d une mèche noire flottante. D une façon générale, les 
réalisations pratiques se sont orientées dans deux voies différentes : 

I ç e ji e des distillateurs de petites dimensions appelés à répondre aux besoins en eau 
potable d’individus isolés ou de petites communautés et même à des besoins plus 
j importants par juxtaposition d’une grande quantité de petits appareils construits en 
série. Celle des distillateurs solaires de plusieurs centaines de m dont la construction 
s’appuie sur des méthodes génie civil. 

0.10- Besoins en énergie pour la distillation de Teau[7] 

L’énergie exigée pour évaporer l'eau est la chaleur latente de vaporisation de l’eau. 
Ceci a une valeur de 2260 kilojoules par kilogramme (kJ/kg). Cela signifie que pour 
produire 1 litre (c’est-à-dire 1kg puisque la densité d'eau est 1 kg/litre) d'eau pure en 
distillant l'eau saumâtre exige un apport de chaleur de 2260 kJ. Cela ne tient pas 
compte de l'efficacité de la méthode chauffante, qui sera moins de 100 /o, ou pour tout 
rétablissement de la chaleur latente qui est rejetée quand la vapeur d'eau est condensée. 

DL11- Conclusion 

Ce chapitre expose un états de l'art des procédés de dessalement de distillation 
thermique et à effet serre et leurs principes de fonctionnement, et leurs performance 

énergétiques. 
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j_i‘r<ormances thermo-optiques du distillateur 


III.1. Introduction 

Ce chapitre reprend 1 essentiel des travaux de modélisation du distillateur à 
multiples effets à effet de serre en se basant sur des bilans énergétiques et des 
propriétés thermodynamiques globaux. Dans cette approche le distillateur est 
décomposé en deux sous distillateurs à simple effet superposés. Chaque sous système 
fonctionne entre deux sources de chaleur chaude et froide caractérisant l'eau saumâtre et 
l'ambiance. Le fonctionnement du distillateur global est régit par une loi linéaire 
exprimant l'efficacité globale en fonction de la température réduite, mesurant le gradient 
de température entre les deux sources de chaleur sur la radiation solaire incidente sur le 
plan d’ouverture permettant la détermination du rendement optique et du facteur des 
pertes globales. Les températures de fonctionnement sont aussi présentées et analysées. 


UI.2. Modélisation des échanges thermiques au niveau du distillateur 
solaire multiple effet de serre 

10.2.1. Description et principe de fonctionnement: 


a- Description 

Le dispositif expérimental (figure. 1IÏ.1) est un distillateur à effet serre à double 
effets. Le système est constitué de deux compartiments superposées lune au dessus de 
l'autre et fonctionnant selon l'effet de serre. Le bac absorbeur de l'unité de base est 
construit en acier inoxydable et peint en noir mat, de lm de surface. Les surf 
latérales des deux unités sont peintes en peinture réfléchissantes pour améliorer la 
réflexion du rayonnement solaire dans les bacs et les saumures. Le bac de base 
20 litres de saumure, et recouvert d'un vitrage ordinaire de 4mm d épaisseur 
supérieure est recouverte en verre ordinaire de 4mm d épaisseur et contien 
d'eau saumâtre. Le système est isolé thermiquement par une couche de polystyrène de 
d'épaisseur. Les niveaux de saumure dans le bac sont fixés par des fl 
rebords du système, sont pratiquées des orifices connectés 
(collecteurs), destinées à alimenter le distillateur en eau saumatre, a nncer les bassins, 
recueillir les distillats et conduire les thermocouples au cours des essais expé 
du distillateur. 

. . rqoutchouc et des mastiques 

Le système est assemblé par des boulons, des jointures 

• > j'oriw soudé éleve de 7U cm 
s dicones. L’ensemble repose sur une structure en comiere 

au-dessus du sol et oriente plein sur site d Adrar. 
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b-Principe de fonctionnement 

Le distillateur solaire à double effets de serre (figure. III.1) est composé de deux 
Mités superposées. Ces unités sont deux sous distillateurs à effet serre. Ainsi, chaque 
sous-systeme est assimile a une boîte noire qui fonctionne entre deux sources de chaleur 
(chaude et froide). Pour simplifier l'étude du système, on peut considérer chaque sous- 
distillateur comme un stege d'une isolation thermique pour le sous distillateur voisin. 



Figure III-l: Distillateur solaire à multiple-effets à effet de serre 

Les principales hypothèses de modélisation sont : 

Chaque sous système est assimilé à une boite noire fonctionnant entre deux sources 
de chaleur (saumure, ambiance). 

les propriétés thermophysiques (transmitivité, absorptivité, conductivité, densité, 
capacité, viscosité,...) des composants de chaque sous système sont indépendantes 
des températures. 

Les facteurs d'échanges par conduction, par convection et par rayonnement sont 
modalises par un seul facteur de pertes thermique global indépendant de la 

température. 

Bilan thermique dans le distillateur [9,10, 111 

La surface d'ouverture du distillateur reçoit une fraction de la radiation solaire Gi. 
Ce « e fraction est transmise par le vitrage et absorbée par la saumure Q* s (source 
cha “de). Une partie de celte énergie est absorbée par l'eau saumâtre et elle est utilisée 
«""me gain utile pour évaporer une quantité d'eau à la température moyenne de la 
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saumure. La quantite de chaleur qui reste est dissipée dans l'atmosphère Q p (source 
froide). L'application du premier principe de la thermodynamique (sortie — entrée 
+$tockage = 0) à un volume de contrôle délimitant les surfaces fonctionnelle du 
système, permet d'écrire l’équation l'équilibre thermique du système avec son 
environnement. 

Au niveau du sous système supérieur: 


Qu.\~Qabs.\ Q P .\ ~ U li A c] {T E j T a ) (III. 1) 


- Au niveau du sous système inférieur 


~Qp,2 ~ ^c,2 a i,2^\ ” U L J A c 2 {r E '2 “ T’a) (IU.2) 


Les fractions de chaleur Q u] , Q u2 représentent respectivement les gains utiles des sous 
distillateurs supérieur et inférieur. 

les quantités de chaleur ,. fi*sont respectivement les flux de chaleur absorbés 
parle saumure dans le sous distillateur supérieur in inférieur. 

Les quantités de chaleur Q p ] , Q p l sont respectivement les flux de chaleur dissipés par 
le sous distillateurs supérieur et inférieur. 


4 p d c 2 sont les surfaces d'ouverture des deux sous systèmes. 

cr (l ,a, 2 représentent respectivement les coefficients d absorption globale 
(rendement optique) de chaque sous système. 

Gp L’irradiation globale sur le plan d'ouverture du sous système supé 

%, U L1 sont respectivement les coefficients des pertes thermiques globales des sous 
distillateurs supérieur et inférieur. 

r„, r„ ,es 'm—. *• —” »“ 

et inférieur. 


L'- La température moyenne de l'ambiance. 
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a* Rendement optique du distillateur 


L'efficacité optique du distillateur supérieur est estimée à partir des fractions de chaleur 
transmise par le vitrage et absorbée par la saumure en additionnant celle qui est 
transmise par le vitrage et par la saumure et absorbée par le vitrage. 


a ;,l - r v,l a E,l + r v,J r £,l ûr v,l (LLL3) 

fl est tout à fait clair que si on néglige le terme * Vt \TE,i a \>,[ ~ 0> on retrouve 
l'absorptance effective d'un capteur plan. 

Par analogie le rendement optique du distillateur inférieur peut être estimé par 
l'expression : 

a ,,2 = ( r v,l - «U l a E,2 + T E,2 a v,2 ) (Hl- 4 ) 

r Tl , r r2 : sont respectivement les transmittances des vitres des deux sous systèmes 
supérieur et inférieur 

a k;L : est le coefficient d'absorption du bac du sous distillateur inférieur. 

r £J ,r f i : représentent respectivement les transmittances des saumures. 

a E j, a E , représentent respectivement les coefficients d'absorption des saumures. 

Tableau III. 1- Paramètres optiques des composantes d'un distillateur solaire simple 

effet[ 1| 


Angle d'incidence du rayonnement 

0-30 

solaire en dégrées 


vitre 

Réflexion % 

5 


Absorption% 

5 


Transmission% 

90 

Nappe 

Réflexion % 

2 

d'eau 

Absorption% 

30 


Transmission% 

68 

Fond 

Réflexion % 

5 

distillateur 

Absorption% 

95 


Transmission% 

0 


Pour un angle d'incidence inférieur à 15 degrés 


a, ,~0.3 et a,, 2 = 0.57. 
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b- Facteur des pertes thermiques 

Le coefficient des pertes thermiques globales est estimé à partir des facteurs des 
pertes thermiques arrières, avants et latérales pour chaque sous distillateur supérieur^ , 

«(inférieur U L 2 similaire au collecteur solaire plan: 


^i.,1 ^ avait,\ ^arriére,\ *U Mt (IU.5) 

= U^ (IU.6) 

Umriirt, i et U itérai.i représente respectivement les coefficients des pertes avant, 

«rière et latéral du sous distillateur supérieur, et 2 , U amitr2 et V laint:i sont 

respectivement les coefficients des pertes avant, arrière et de latéral du sous distillateur 
inférieur. 

U facteur de perte global sera alors: 

Vl=u u +u L ' 2 (in.7) 

IDiJ. Le rendement thermodynamique du distillateur 

Le rapport entre le gain utile et l’énergie solaire incidente sur le plan d'ouverture du 
Milateur durant une de temps donne l'efficacité thermique globale des sous-systèmes 
comme suit : 



(IH.8) 

(ffi.9) 


fyi'fy .2 représente respectivement les efficacités thermiques des sous-systèmes 
supérieur et inférieur. 

m 2 sont respectivement les débits d'eau distillée dans les sous-systèmes supérieur et 
inférieur. A c i , A c 2 : sont les surfaces pertinentes. E g ] : l'énergie solaire incidente . 

L y} en kj/kg: sont respectivement chaleur latente de vaporisation de l’eau dans les 
Gu$>$ystèmes supérieur et inférieur. 
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/. vl =4.1 «[883-0.668(7', ,+273.15)] (UI.10) 
/., 2 = 4.18(883-0.668(7', ;2 + 273.15)] (UI.ll) 


La combinaison des équations (III. 1), (III.2) avec les équations (III.8) et (III.9), donne 
l'expression des efficacités linéaires des sous systèmes sous les formes: 


u i,i 


tfd.l — a i.2 U L. 


£.1 1 a 

\ J 
V G, ) 


III 

1 

to 

(ÏU.12) 

sA 2 - B 2 T r2 

(111.13) 


les constantes Ai, Bj, A 2 et B 2 sont déterminés expérimentalement. 


T 0 : est la température réduite des sous distillateur supérieur. 
ï ti : est la température réduite des sous distillateur inférieur. 


G,: l'éclairement global. 


Selon le principe de superposition, le rendement global du distillateur à multiples effet 
est déduit à partir des rendement des deux distillateur sous la forme: 


Vd ~ — ~pr Â - ~( a t, 1 +a t, 2 )“ 

UA,. 


U 


L, 1 


(7'£,l- 7 d.„ (Te,2<) 

-*■ VI J - ~z 


(in.i4) 


les gradients T t , - T a ' T h2 - T a sont pratiquement constantes ce qui permet de déduire 
la loi de fonctionnement du distillateur à multiple-effets: 

=K, +«fcMÏ4 > i +64^^ 011.15) 

G[A 4 4 
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111.2.4. Mode opératoire 

Cette étude à été effectuée par l'Unité de Recherche en Lnergie Renouvelable en 
Milieu Saharien d'Adrar (latitude 27.88 °N. longitude = 0.28° Ouest et une altitude de 
264 m ). Le distillateur solaire est exposé au soleil aux alentours de midi temps solaire 
vrai pour un angle d'incidence <15 degrés, une vitesse de vent <2 m/s, un éclairement 
de 1000 W/m 2 et température ambiante <30. Le rayonnement solaire incident sur la 
surface d'ouverture du distillateur solaire a été mesuré par un pyranomètre Kipp & 
Zonnen de seconde classe. La vitesse du vent a été mesurée avec un anémomètre de 
précision. Les températures ambiantes et des saumures sont mesurées par des 
(homocouples de type K reliés à un dataloger (enregistreur de données) de 21 voies de 
type (Fluke 2635A) à mémoire de stockage incorporé et les données sont stockées 
toutes les 5 minutes pendant les essais. Les quantités d’eau distillée sont mesurées par 
deux vases grad ués ( fi g u re 111.2). 



(D B*ia 

(D BHm 

® GImcmw 
0 IimMm 
® Fl* 

0 S*pJy«tfk« 
0 MNOilIct 
® CdkâMtta 



Figure IH-2: Schéma du dispositif expérimental 
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Performances thermo-optiques du distillateur 


IIL3- Résultats et discussions 


«000 -1 



100 woouoo wooiooo noo»«aooï*aojioejo«o»ooMoo 
Tafnpaloca<(h) 



ioo o.co moi 2.001 « oow.ooiiooM.ooa ow.ooai mtjooo n ooaoo 
Twnÿtoca(h) 


(a) 


(b) 


Figure III-3: (a) Variation des températures de vitrage des deux sous distillateurs, 
(b)variation temporelle du gradient de température de vitre (Tv inr Tv Mp ) 


La courbe (a) représente la variation temporelle de température des vitres dans les 
deux sous distillateurs inférieur et supérieur. On remarque une élévation des 
températures des deux vitres en fonction de l’augmentation de l’éclairement solaire 
incident, la réduction des déperditions thermiques. Cette élévation devient maximale 
vers 13 à 14 heures temps local, avec un maximum de 55°C pour la vitre de sous 
distillateurs inférieur. En même temps, la température de la vitre du sous système 
supérieur atteint un maximum de 44°C. Elle est inférieur à celle de sous système 
inferieur à cause de pertes de chaleur par convection vers l’ambiance sous l’action du 
vent. L'après midi TSV, on remarque une diminution des températures des vitres de 
chaque sous systèmes due à la diminution de l'éclairement et l’augmentation des pertes 
thermiques vers l'ambiance jusqu’à l’équilibre total du système avec son environnement. 
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(*) (b) 

Figure ni-4: (a) Variation du température de saumure dans les deux sous distillateurs 
inférieur et supérieur, (b): Variation temporelle de l’écart température de saumure (T Ei-f 


La courbe (a): représente la variation des températures des eaux saumâtres des sous 
systèmes inférieur est supérieur. Le gradient de température entre ces deux sous 
systèmes augmente avec l'augmentation du rayonnement solaire pour dépasser les 20°C 
en raison de l’isolation thermique avant et arrière du distillateur de base. Ainsi, la 
production se trouve amélioré. Toutefois les fluctuations sont similaires à celles des 
vitrages et l'interprétation est donc semblable. 



(a) (b) 

Figure 111-5: (a): Variation des températures des bacs des deux sous distillateurs inférieur 
et supérieur, b: Variation temporelle de l'écart température des bacs (Tbi„r — Tb SU p) 
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(a) Sous distillateur inférieur 


(b) Sous distillateur supérieur 



(c) Les deux sous distillateurs 

Figure III-6: Température de fonctionnement des sous distillateurs 

La figure III-6 montre la variation de température de chaque composante de sous 
distillateur avec sont environnement en fonction du temps. Les niveaux de température 
sont conditionées par le rayonnement solaire, les propricts thermophysiques et les pertes 
thermiques et optiques. Les températures des bacs sont les plus éleveés et celles ds 
vitrages sont les plus faibles. Les courbes présentent des augmentations progressives de 
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Pci tormaïuTs tlnrmo-optiques du distillateur 


températures en fonction du rayonnment solaire du lever de soleil à 12 tsv, ensuite une 
relative constance durant 2 à 3 heure indiquant une variation nulle de l'énergie interne 
du système, c’est une phase d'équilibre entre les pertes thermiques et les gains utiles 
suivi d'une décroissance progressive indiquant la phase de déchage thermique 
augmentaion progressive des pertes thermique le soir pour arriver à un équilibre 
thermodynamique du système avec le milieu extérieur en fin de journée. 



(a) Sous distillateur iuférieur 


(b) Sous distillateur supérieur 


Figure I1I-7: Variation des débits de distillât dans les sous distillateurs en fonction du 

gradient (T E -T,) 



(b) Sous distillateur inférieur (b) Sous distillateur supérieur 


Figure I1I-8: Variation des débits de distillât en fonction de l’éclairement incident 
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Les débits des distiilats dans les deux sous distillateurs augmentent avec 
l'augmentation des gradients des températures entre les sources chaudes (saumures) et 
les sources froides (ambiances) (figure II 1-7). Également les débits augmentent avec 
l'augmentation de la radiation solaire incidente (figure III-8). Selon les droites de 
corrélation, l'augmentation dans le sous distillateur inférieur est plus importante que 
celle dans le sous système supérieur pour les raison précitées. 



(a) Sous distillateur inférieur (b) Sous distillateur supérieur 

Figure III-9: Caractéristiques thermo-physiques des sous distillateurs 

Les courbes de la figure III-9, représentent les caractéristiques principales des sous 
systèmes; exprimant le rendement de chaque sous distillateur en fonction de la 
température réduite. Le nombre de point de mesure choisis répond relativement aux 
conditions de caractérisation des capteurs plans la plus proche du standard ASHRAE 
(1000 W/m 2 , 25°C 0 à 3 m/s, incidence inférieure à 15 degrés). L'approximation du 
nuage de points par la méthode des moindes carrées donne une droite. L'intersection de 
la droite avec l'axe des abscisses dorme les facteurs des pertes thermiques globales qui 
(UL,i~3.6W/m 2 C et Uj , t 2 ~ 2,6 w/m 2 C) et l'intersection avec l'axe des coordonnées donne 
les rendements optiques qui sont (a,,i~24% et a t ,2~56%). Le facteur des pertes 
thermiques du distillateur supérieur est relativement réduit car le sous système 
compense les pertes thermiques par la chaleur de condensation et de rayonnement du 
sous système inférieur. 
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La comparaison des résultats obtenus avec les valeurs estimées par le calcul direct à 
partir des propriétés thermo-physiques des sous systèmes donne (Ul,i = 3,5 W/m C et 
a t ,i~30%) et (l)| 2 2 W/rrf(' et a,,2-57%) avec des erreurs respectives de (2.78% et 25 

%) pour le système inférieur et de (23.01% et 1.78 %) pour le système supérieur. 
Cependant le distillateur à multiple effets est caractérisé par un rendement optique de 
87% et un facteur de pertes de 5,5 W/m 2o C. En tenant compte des hypothèes 
simplificatrices, des ererurs accidentelles et systématiques ces résulats sont relativement 
accepatbles. 

ffl.4- Conclusion 

Ce chapitre a expose une modélisation simplifiée de distillateur à effet de serre 
doubles effet aussi il a permet de déterminer les caractéristiques principaux de 
système. Les résultats obtenus sont illustrés dans le tableau suivant: 

Tableau III.2. Les performances de système 


Système 

Absorbité 

effective 

[%] 

Absorbité 

Mésuré 

(%] 

■ 

■ 

Coefficient 

de perte 

effectif 

[W/m 2o C] 

Coefficient 

de perte 

mesuré 

[W/m 2 °C] 

Erreur 

[%] 

Sous distillateur 

supérieur 

30 

24 

20 

3.5 

3.6 

2.78 

Sous distillateur 

inférieur 

57 

56 

1.75 

2 

2.6 

23.01 

Distillateur multiples 

effet 

87 

80 

Q 

5.5 

6.2 

25.79 


Le rendement et le facteur des pertes thermiques du sous distillateur supérieur est 
relativement réduit car le sous système compense les pertes thermiques par la chaleur de 
condensation et de rayonnement du sous système inférieur. 
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Conclusion generale 


Conclusion générale 

Cette modeste approche de modéli salon et de caractérisation d'un distillateur thermo sol aire à 
multiple effets a permis de mieux comprendre les sources d'énergie solaire les bilans 
énergétiques mis en jeu dans le système avec l'ambiance et les propriétés thermodynamiques 
d'évaluation. Elle a permis d’appliquer une méthode globale simple basé sur des bilans 
énergétqiues et des propriéts thermo-optiques pour caractériser le système en se basant sur un 
minimum de grandeurs physiques et de moyens de mesure. Cette procédure semi-empirique a 
permis de donner les deux pramètres clés de design du système : le rendement optique et le 
facteur des pertes thermqiue globales dans des conditions les plus proches des standards 
internationales des tests des capteurs solaires plans. Les résulatts obtenus sont relativement 
acceptabels en tenant compte des hypothèses simplificatrices. À surface égale le distillateur 
thermosolaire à multiples effets permet d'améliorer la production de distillât de 10 à 24% par 
rapport au distillateur à simple effet. Cette approche de caractéridsation constitue une base pour 
la modélisation fine des distillateurs à effet serre à simple et à multiples effet à grande échelle, 
elle aussi cirucial pour le dimensionnment et l’optimisation des prototypes expérimenatux. 

Perspectives 

- Amélioration de la qualité de caractérisation en tenant compte des fluctuations des paramètres 

thermo-optiques en fonction de la température. 

- Introduction de la notion de constante de temps pour l'optimisation de l'inertie thermique du 

système. 

- Proposition de nouvelles configurations de distillateurs thermosoalires en sépara» 

l'évaporateur du condenseur. 
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